Coherent radial breathing like phonons in graphene nanoribbons 
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We have developed a microscopic theory for the generation and detection of coherent phonons in 
armchair and zigzag graphene nanoribbons using an extended tight-binding model for the electronic 
states and a valence force field model for the phonons. The coherent phonon amplitudes satisfy 
a driven oscillator equation with the driving term depending on photoexcited carrier density. We 
examine the coherent phonon radial breathing like mode amplitudes as a function of excitation 
energies and nanoribbon types. For photoexcitation near the optical absorption edge the coherent 
phonon driving term for the radial breathing like mode is much larger for zigzag nanoribbons where 
transitions between localized edge states provide the dominant contribution to the coherent phonon 
driving term. Using an effective mass theory, we explain how the armchair nanoribbon width changes 
in response to laser excitation. 

PACS numbers: 63.22.-m, 73.22.-f, 78.67.-n 
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I. INTRODUCTION 

Excited state lattice vibrations in carbon nanotubes 
have been studied with coherent phonon (CP) spec- 
troscopy^— Using CP spectroscopy with pulse shap- 
ing techniques, radial breathing mode (RBM) coher- 
ent phonons in chirality-specific semiconducting single- 
walled carbon nanotubes have been studied experimen- 
tally^— The CP signals are resonantly enhanced when 
the pump photon energy coincides with an exciton 
resonance, and provides information on the chirality- 
dependence of light absorption, phonon generation, and 
phonon-induced band structure modulation. 

Recently we developed a microscopic theory for the 
generation and detection of coherent phonons in single- 
walled carbon nanotubes in CP spectroscopy experi- 
ments^— We used the Heisenberg equation to obtain a 
driven oscillator equation for the coherent phonon am- 
plitudes and found that the driving function depends ex- 
plicitly on the time-dependent photoexcited carrier dis- 
tribution functions. Comparing theory and experiment 
we find that our model predicts overall trends in the rel- 
ative strength of the RBM coherent phonon signal both 
within and between different tube types. 

In a followup paper we studied the chirality depen- 
dence of coherent phonon amplitudes in single-wall car- 
bon nanotubes over a large range of chiralities using an 
effective mass approach^ By examining the ^-dependent 
electron-phonon interaction in the effective mass approx- 
imation, we were able to explain why some nanotubes 
start their coherent RBM diameter oscillations by ini- 
tially expanding while others start their RBM diameter 
oscillations by initially shrinking their diameters, a fact 
recently observed experimentally^ In many solids, lat- 
tices tend to expand when photoexcited by ultrafast laser 
pulses in accordance with the Franck- Condon principle. 
As pointed out in Ref. this is not the case for RBM 
diameter oscillations in carbon nanotubes where the di- 
ameter can either expand or contract depending on the 



nanotube chirality and photoexcitation energy. 

The electronics industry is exploring device technolo- 
gies based on carbon nanotubes, graphene, and graphene 
nanoribbons (GNRs). 10 Field-effect transistors based on 
GNRs have been demonstrated 11 and it is now possible 
to fabricate GNRs with atomically precise widths using 
a number of methods*^— An understanding of the elec- 
tronic and transport properties of GNRs is essential in 
realizing device applications for GNRs. 18 In particular it 
is important for characterization and transport modeling 
to have a good understanding of the electron-phonon in- 
teraction and lattice vibrations in GNRs. It is hoped that 
coherent phonon spectroscopy will prove useful in charac- 
terizing graphene and graphene nanoribbons in addition 
to carbon nanotubes. 

In this paper, we extend our coherent phonon the- 
ory to the case of unpassivated zigzag and armchair 
nanoribbons. In our model, we calculate electronic states 
for the 7r electrons in an extended tight binding model 
(ETB)^ii^ while the phonon modes are treated in a va- 
lence force field framework] 20 ! 21 Details concerning the 
computation of electronic states and phonon modes can 
be found in Appendices [A] and [B] respectively. In our 
formalism, we incorporate the electron-phonon interac- 
tion^ 2 - the optical matrix elements^ and the interaction 
of carriers with a classical ultrafast laser pulsed For sim- 
plicity, we neglect the many-body Coulomb interaction. 
In carbon nanotubes, we found that there were exactly 
four CP active phonon modes independent of the tube's 
chirality. The most easily observed CP active mode in 
nanotubes is the RBM, which has the lowest frequency. 
In zigzag and armchair nanoribbons, on the other hand, 
the rotational degree of freedom about the nanotube axis 
is lost and the number of CP active phonon modes in 
zigzag or armchair nanoribbons is equal to the number 
of AB carbon dimers in the nanoribbon translational unit 
cell. The most easily observed of these CP modes is the 
one with the lowest frequency, namely the radial breath- 
ing like mode (RBLM). 
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FIG. 1: Schematic diagram showing lattice structures and 
translational unit cells for armchair (aGNR) and zigzag 
(zGNR) graphene nanoribbons. The width of the nanorib- 
bons is W. 



II. THEORY 

We consider two types of nanoribbon, namely zigzag 
and armchair ribbonsi 24 ' 25 The lattice structure for arm- 
chair and zigzag graphene nanoribbons is shown schemat- 
ically in Fig. [TJ These ribbons are denoted N a b zGNR 
and N a b aGNR respectively where N a b is the number of 
AB carbon dimers in the translational unit cell. In zigzag 
ribbons, the length L of the translational unit cell is a 
and the width W of the ribbon is (N a b — 1) where 
a = 2.49 A is the hexagonal lattice constant in graphene. 
In armchair ribbons, the translational unit cell length is 
\fia and the ribbon width is (N a b — 1) \a>- Note that in 
zigzag and armchair ribbons with the same number of 
atoms per unit cell, the area of the unit cells are equal. 
Note that we do not allow for the relaxation of the bond 
lengths at the surfaces. 



A. Coherent phonon generation 

In coherent phonon spectroscopy, the coherent phonon 
modes that are usually excited are the modes with 
wavevector q = whose amplitudes satisfy a driven os- 
cillator equation 



d 2 Q m (t) 
dt 2 



U 2 m Q m {t) = S m {t), 



(1) 



where m denotes the phonon mode and u> m = uj m (q = 0) 
is the frequency of phonon mode to at q = 0. There 
is no damping term in Eq. (fTJ) since anharmonic terms 
in the electron-phonon Hamiltonian are neglected. We 
solve the driven oscillator equation subject to the initial 
conditions Q m (0) = Q m (0) = 0. The driving function 
S m (t) is given by 

Sm (t) = -^J2 M n(k) (/»(M) - fn(k)) ■ (2) 
n k 

where f n (k,t) arc the time-dependent electron distribu- 
tion functions and are the initial equilibrium elec- 
tron distribution functions. Here n labels the electronic 



state and k is the electron wavevector. The electron- 
phonon matrix element M™(fc) = M%g nk where M™ q k ,. nk 
is the deformation potential electron-phonon matrix ele- 
ment in our ETB model. Details concerning calculation 
of this matrix element can be found in Appendix [Cl 

The driving function S m (t) depends on the photoex- 
cited electron distribution functions which can be calcu- 
lated in the Boltzmann equation formalism taking pho- 
togeneration and relaxation effects into account. In CP 
spectroscopy, an ultrafast laser pulse generates electron- 
hole pairs on a time scale short in comparison with the 
coherent phonon period. The observed CP signal is pro- 
portional to the power spectrum of Q m (t). For simplicity 
we ignore relaxation effects and retain only the rapidly 
varying photogeneration term in the Boltzmann equa- 
tion. This model works best when the carrier relaxation 
time is greater than the coherent phonon period. The 
ultrafast pump pulse gives an impulsive kick to Q m (t) 
setting the coherent phonon oscillations in motion. The 
subsequent slow relaxation of S m (t) shifts the coordinate 
about which Q m (t) oscillates but the coherent phonon 
oscillator is able to follow this shift adiabatically with 
negligible change in the oscillation amplitude and the cor- 
responding power spectrum at the oscillation frequency. 
Thus neglecting carrier relaxation has a negligible effect 
on the computed CP signal. 

The photogeneration rate is obtained from Fermi's 
golden rule and the equation of motion for the distri- 
bution functions on time scales short in comparison to 
the relaxation time is 



df n (k) = 8^ 2 e 2 u(t) 

dt h v? (hbj) 2 



m 



x(/n'(M)-/»(M)) 5(AE nn ,(k)-?kj), (3) 

where AE nn >(k) — \E n (k) — E n /(k)\ are the k dependent 
transition energies, hu is the pump energy, u(t) is the 
time-dependent energy density of the pump pulse, e is 
the electron charge, mo is the free electron mass, and 
n g is the index of refraction in the surrounding medium. 
The pump energy density u{t) is related to the fluence 
by F — J dt u(t) (c/rig). The pump energy density is 
taken to be a Gaussian with an intensity full width at 
half maximum (FWHM) of r p which we define as the 
pump duration. The optical dipole matrix element is 
given by 



Pnn' (k) — 



y/2m 



dv < fc (r) W„, fe (r), (4) 



where e is the complex electric polarization vector of unit 
length and ip n k (r) are the n electron tight-binding wave- 
functions defined in Appendix [3] With the aid of the n 
orbital expansions of the electron wavefunctions we can 
evaluate the optical dipole matrix element in Eq. (jj) ana- 
lytically. To account for spectral broadening of the laser 
pulses we replace the delta function in Eq. ([3]) with a 
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Lorentzian lineshap> 



,20 



5{AE - fvjj) 



r p /(27r) 



(AE - hu) 2 + (r p /2) 2 



(5) 



where T p is the FWHM spectral linewidth of the pump 
pulse. 

In addition to the photogeneration rate, there are also 
carrier relaxation effects to consider. In Ref. [27] Dawlaty 
et al. measured carrier relaxation times in photoexcited 
graphene using pump-probe spectroscopy and found two 
relaxation time scales for relaxation of photogenerated 
carriers. An initial fast relaxation transient with a re- 
laxation time ranging from 70 to 120 fs is followed by 
a slower relaxation with a relaxation time ranging from 
400 to 1700 fs. These effects can be modeled in a phe- 
nomenological relaxation time approximation which adds 
a term to the photogeneration rate in Eq. (j3|) of the form 



d.Uk) 
dt 



relax 



f n (k,t) - f°(k) 



(6) 



where t t is the phenomenological relaxation time and 
fn(k) are the initial carrier distribution functions in ther- 
mal equilibrium. 

In our tight-binding model we only include 7r bands. 
The a bands in graphene have a direct gap at the T point 
of around 6 eV.— In nanoribbons the a band direct gap 
will be even larger due to quantum confinement effects. 
The laser pulses discussed in this work are restricted to 
photon energies fku < 5 eV so we need not consider co- 
herent phonon generation due to photoexcited carriers 
from the a bands. 

From the coherent phonon amplitudes, the time- 
dependent macroscopic displacements of each carbon 
atom in the nanoribbon are given by 



(7) 



The imaginary part of the nanoribbon dielectric func- 
tion is obtained from Fermi's golden rule 



E 



dk 



WL z (Huj) 2 \m 

(/n(*0 - /n'(fc)) - E n (k) - hu 



Pnn' X 



(9) 



where W is the nanoribbon width and L z — 3.4 A is 
the nanoribbon thickness taken to be the interlayer dis- 
tance in graphite. We replace the delta function in 
Eq. (0) with a broadened Lorentzian spectral lineshape 
with a FWHM of T s . The distribution function f n (k) 
and bandstructure E n (k) are time-dependent. The time- 
dependence of f n {k) comes from the photogeneration of 
carriers described by the Boltzmann carrier dynamics 
and the time-dependence of E n (k) comes from variations 
in the carbon-carbon bond lengths due to the macro- 
scopic coherent phonon induced atomic displacements 
in Eq. ([7]). This time-dependent deformation of the 
nanoribbon bond lengths alters the tight-binding Hamil- 
tonian and overlap matrix elements in the extended tight- 
binding model.- Note that to first order in the lattice dis- 
placements the energies E n (k) vary with time while the 
tight-binding wavefunctions and optical matrix elements 
Pnn'(k) do not. 

In coherent phonon spectroscopy, excitation of coher- 
ent phonons by the pump modulates the optical prop- 
erties of the nanoribbons giving rise to a transient dif- 
ferential transmission signal. In our model we take the 
theoretical CP signal to be proportional to the power 
spectrum of the transient differential transmission after 
background subtraction. We compute the power spec- 
trum using the Lomb periodogram algorithm described 
in Ref. M. 



III. ZIGZAG NANORIBBON RESULTS 



A. Bandstructure and absorption spectra 



where e™ = e" l Jq = 0) is defined in Appendix IB"1 



B. Coherent phonon detection 

In coherent phonon spectroscopy a probe pulse is used 
to measure the time- varying absorption coefficient of the 
nanoribbon. The time-dependent absorption coefficient 
is given b y 26 ' 28 



a{hui, t) 



Tig HC 



(8) 



where E2(hui,t) is the imaginary part of the time- 
dependent dielectric function evaluated at the probe pho- 
ton energy Hoj. 



In zigzag nanoribbons, the nanoribbon translational 
unit cell length is L = a = 2.49 A and the nanoribbon 
width is W = (N a b — 1) ^a. We consider in detail a 
zigzag nanoribbon with 7 dimers per translational unit 
cell (7 zGNR). In this case N a b = 7 and the width of the 
ribbon is W — 12.94 A. The presence of the edges in the 
nanoribbon qualitatively alters the electronic and optical 
properties. 

The computed bandstructure E n {k) and density of 
electronic states DOS(E) for 7r bands in 7 zGNR nanorib- 
bons is shown in Fig. [Ha) and (b) . The Brillouin zone is 
one dimensional with |fc| < n/L. There are total of 2N a b 
or fourteen bands. The seven bands with negative en- 
ergy (E n (k) < 0) are the valence bands and the seven 
bands with positive energy are the conduction bands. 
From Fig. [2fa) we note that the conduction and valence 
bands are asymmetric about E — 0. This comes from the 
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FIG. 2: (color online) (a) Computed ETB electronic bands 
for 7T electrons in a zigzag nanoribbon with 7 atoms per trans- 
lational unit cell, (b) Corresponding electronic density of 
states. The density of states per carbon atom for the valence 
bands are shaded in yellow. 



atomic overlap matrix elements in the ETB formalism 
which is not present in the simple effective mass model. 
We see a strong degeneracy in the bands at k = ± ir/L 
and note that there are two partially degenerate bands 
near E = 0. These bands correspond to localized edge 
states for \k\ > 2ir/3L where k = ± 2ir/3L are the Dirac 
points K and K' . For |fc| < 2ir/3L the states penetrate 
into the interior of the ribbon as |fc| — > 0. These edge 
states are peculiar to zigzag nanoribbons and have been 
studied using effective mass, tight-binding and ab ini- 
tio theories^^i^. The remaining bands are delocalized 
zone folded quantum confined bands the lowest of which 
have parabolic minima near the Dirac points K and K' 
located at k = ± 2ir/3L. The density of states shown in 
Fig. HJb) has a sharp peak near E — due to the disper- 
sionless localized edge bands in the vicinity of A: = ±ir/L. 
As we move away from E = the density of states con- 
tains a series of von Hove singularities at the extrema 
of the electronic bands. As N a i, increases, the density of 
states approaches DOS(E) in planar graphene. For large 
N a b the edge state contribution to DOS(E) decreases as 
l/N ab M 

The absorption spectrum is shown in Fig. [3ja) for the 
7 zGNR nanoribbon. For now we assume the light is lin- 
early polarized with the electric polarization vector par- 
allel to the nanoribbon's length. Later we will study the 
in-plane polarization dependence. The nanoribbon is as- 
sumed to be undoped at a temperature of 300 K. The 
Fermi level E-p at 70.16 meV is indicated by the black 
dashed line in Fig. [3Jb). The absorption coefficient is 
computed using Fermi's golden rule (Eq. (|8]l) assuming 
Lorentzian lineshapes with a FWHM of 0.2 eV. The to- 
tal absorption spectrum is shown in Fig. [21a) as a thick 
black line. This is a sum of contributions from several 



FIG. 3: (color online) (a) Band edge absorption for 7 zGNR 
nanoribbons with electric polarization vector parallel to the 
ribbon. Total absorption is the sum of absorption due to sev- 
eral transitions, (b) Band diagram showing the Fermi level 
Ef and transitions involved in the three lowest absorption 
peaks. Blue lines are valence bands and red lines are conduc- 
tion bands. Solid lines (n even) are even parity states and the 
dash-dotted lines (n odd) are odd parity states. 



transitions. The lowest lying absorption peak at 1.65 eV, 
labeled V2C1, is due to transitions between the low lying 
second hole band V2 and the higher lying first electron 
band CI. This transition is indicated by a vertical arrow 
in Fig. EJb) whose length is the absorption peak tran- 
sition energy, i.e. 1.65 eV. From Fig. EJb), we see that 
the V2C1 peak comes from transitions between states 
near the K and K' points in the Brillouin zone. The 
initial states are quantum confined hole states V2 and 
the final states are the localized electron edge states CI. 
A second broad absorption transition VI C2 peaking at 
2 eV comes from transitions between the localized hole 
edge states VI and the second quantum confined electron 
states C2. There is a strong peak in the absorption spec- 
trum at 2.23 eV due to the intraband transition V3V1 
between V3 quantum confined states near the K and K' 
points and the localized VI hole edge states. 

The absorption spectrum in the 7 zGNR nanoribbon is 
qualitatively different from the absorption spectrum seen 
in nanotubes. A mirror plane runs down the center of the 
nanoribbon so the electronic states are either symmetric 
or antisymmetric about the mirror plane. The bands 
with even parity in Fig. |3jb) are solid lines and the odd 
parity bands are shown as dash-dotted lines. The allowed 
transitions for light polarized parallel to the nanoribbon 
axis are between bands with the same parity. In partic- 
ular the interband An = selection rule which holds for 
nanotubes docs not hold for nanoribbons. These conclu- 
sions are in agreement with results obtained within the 
effective mass formalism^ The edge states play an im- 
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FIG. 4: (color online) (a) Phonon dispersion relations for 
a 7 zGNR nanoribbon calculated in the valence force field 
model. The out-of-plane modes are red dash-dotted lines 
while the remainder are in-plane modes. The band associ- 
ated with the q = RBLM mode is shown as a thick blue 
line, (b) Density of modes per carbon atom for the dispersion 
relations in (a). 



portant role in the band edge transitions and the absorp- 
tion spectrum is sensitive to the position of the Fermi 
level. In Fig. [3^b) we see that the VI edge states in 
the vicinity of the Dirac points lie above the Fermi en- 
ergy. Since these states are empty, electrons in the lower 
bands can be photoexcited into these VI states and this 
accounts for the strong intraband V3V1 absorption peak 
seen at 2.23 eV. 

Optical selection rules and absorption in zigzag 
nanoribbons have been studied by other authors using 
nearest-neighbor tight-binding approaches . 32 ' 33 They ne- 
glect overlap matrix elements so that the conduction and 
valence bands are symmetric about E = 0. For symmet- 
ric valence and conduction bands the absorption peaks 
V1C2 and V2C1 are degenerate which disagrees with 
what we find taking overlap matrix elements into ac- 
count. Our edge state bands VI and CI are not com- 
pletely dispersionless and in undoped 7 zGNR nanorib- 
bons we find E-p lies below the VI band for a range 
of k values in the Brillouin zone. This gives rise to a 
strong V3V1 absorption feature which is not predicted 
in Refs. Eland EE 



B. Phonon dispersion relations 

Phonon dispersion relations for a 7 zGNR nanoribbon 
obtained in our valence force field (VFF) model is shown 
in Fig. Eta). The dispersion relations and normal mode 
displacements are obtained by solving the dynamical ma- 
trix eigenvalue problem.— Since the force constant tensor 
is block diagonal, the in-plane and out-of-plane modes 



can be treated independently. There are 2 N a b = 14 out- 
of-plane modes and these are shown as red dash-dotted 
lines. We find that for normally incident pump pulses, 
none of these out-of-plane modes can be coherently ex- 
cited and so they are not of interest to us. The remaining 
4 Nat, = 28 modes are in-plane modes and are shown as 
solid lines. We find that N a }, = 7 of the in-plane q = 
modes can be coherently excited. We focus our atten- 
tion on the lowest lying of these modes, the so-called 
radial breathing like mode (RBLM). The RBLM mode is 
a stretching mode in which the width of the nanoribbon 
breathes in and out and is somewhat analogous to the 
radial breathing mode (RBM) seen in nanotubes. The 
band associated with the RBLM is shown in Fig. HJa) as 
a thick blue line. The phonon energy of the RBLM coher- 
ent phonons is Hlu = 29 meV which corresponds to an os- 
cillation period of 0.141 ps. The density of phonon modes 
is shown in Fig. EJb). Apart from numerous von- Hove 
singularities, the density of phonon states in Fig.Etb) re- 
sembles the density of phonon states in bulk graphene. 21 
The RBLM phonon energies in zGNR nanoribbons ob- 
tained from our VFF model can be fit to fvu) = AW P + B 
where W is the zigzag nanoribbon width in Angstroms 
and haj is the q = RBLM phonon energy in eV. The 
fitting parameters A = 0.25247, P = -0.82773 and 
B = —0.00136 are obtained from a least squares fit to our 
VFF results for 3 < N a b < 15 corresponding to widths 
in the range 4 A < W < 30 A. 



C. Generation of coherent phonons 

In a typical simulation, we excite coherent RBLM 
phonons with a 20 fs Gaussian laser pulse. In our sim- 
ulations the fluence F is taken to be 5 x 10 4 J/cm 2 and 
the pump pulse is linearly polarized with the electric po- 
larization vector parallel to the nanoribbon length. The 
pump photon energy is 2.24 eV which coincides to the 
V3V1 absorption peak shown in Fig. [Ha). The pump 
energy density u(t) is plotted against the left axis of 
Fig. [5jb) and the photoexcited carrier density per unit 
length is plotted against the right axis. The time de- 
pendent distribution functions and photogenerated car- 
rier density are obtained by integrating the photogen- 
eration rates in Eq. (J3|) using the equilibrium distribu- 
tion functions as initial conditions. In integrating this 
equation, we assume a FWHM pump spectral width of 
r p = 150 meV which appears in the broadened delta 
function. Using the time-dependent distribution func- 
tions, the RBLM coherent phonon driving function S(t) 
can be obtained from Eq. ([2J. The driving function 
in our example is shown in Fig. [5ja). It initially rises 
sharply during the pump phase and reaches a peak value 
S'max = 94.6 ps~ 2 . Given S(t) we can solve Eq. ([1]) to 
obtain the RBLM coherent phonon amplitude Q(t). Pos- 
itive Q(t) for the RBLM mode corresponds to an expan- 
sion of the nanoribbon width while negative Q(t) corre- 
sponds to a contraction. For coherent RBLM phonon os- 
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FIG. 5: (color online) For 7 zGNR nanoribbon excited by 
Gaussian laser pulse with polarization vector parallel to rib- 
bon length we plot (a) coherent RBLM phonon driving func- 
tion S(t) and (b) pump photon energy density (left axis) and 
photogenerated carrier density (right axis). For the assumed 
pump pulse, the photon energy u, fluence F, pulse duration 
t p , and FWHM spectral linewidth Y v are indicated in (b). 
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FIG. 6: (color online) For 7 zGNR nanoribbon excited by 
a 20 fs Gaussian laser pulse with polarization vector parallel 
to ribbon length we plot (b) coherent RBLM phonon driv- 
ing function S(t) for several relaxation times. Solid curves 
are labeled with the relaxation time measured in picoseconds. 
The dashed blue curve is the relaxationless limit. In (a) the 
dashed blue curve shows the coherent phonon amplitude Q(t) 
in the relaxationless limit and the solid red curve shows Q(t) 
assuming a 10 fs carrier relaxation time. 



cillations the nanoribbon width initially expands or con- 
tracts depending on whether S max is positive or negative. 

Carrier scattering effects can also effect the strength of 
the coherent phonon oscillations especially if the scatter- 
ing times are short compared to the oscillation period. 
A detailed treatment of carrier scattering is beyond the 
scope of this article. However, we note that carrier scat- 
tering by itself does not necessarily lead to a change in 
the driving function S m (t). According to Eq. [2 scatter- 
ing affects S m (t) only if the carrier scatters to a state 
that has a substantially different carrier-phonon matrix 
element. For instance, in GaAs, an electron scatter- 
ing within the conduction band (either via incoherent 
phonons or electron-electron scattering) does not change 
the driving function since the electron-phonon matrix ele- 
ments are similar within the conduction band. We would 
expect that in graphene nanoribbons, scattering within a 
given band would have less of an effect than scattering 
between different bands. 

We can estimate the effects of carrier dynamics on 
these results using a simple relaxation time approxima- 



tion model for recombination of photogenerated carriers. 
We considered relaxation times ranging from 10 fs to 2 ps 
and our results are shown in Fig. [6l After an initial sharp 
rise due to photoexcitation by the pump, S(t) decays ex- 
ponentially at long times due to carrier relaxation as seen 
in Fig. [6(b). We then solved the driven oscillator equa- 
tion (Eq. [TJ for the coherent phonon amplitude Q{t). In 
Fig-EJa) we plot Q{t) assuming a relaxation times of 10 fs 
(solid red curve) and in the relaxationless limit (dashed 
blue curve). 

We note that there is a shift in phase of about 7r/2 be- 
tween the two cases. For the 10 fs carrier relaxation time, 
we get a sinusoidal oscillation, while for the no relaxation 
case, we find that Q(t) cx 1 — cos(u)t). These are the lim- 
iting cases for a driving function that is proportional to 
a delta function and a step function respectively. This is 
discussed in more detail in Ref. [34]. 

The CP signal is proportional to the square of the am- 
plitude of the Q(t) oscillations, so to estimate the effect 
carrier relaxation has on the computed RBLM CP sig- 
nal, we solved the driven oscillator equation and found 
the amplitude of Q(t) at long times as a function of the 
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FIG. 7: (color online) For 7 zGNR nanoribbon excited by 
Gaussian laser pulse with polarization vector parallel to rib- 
bon length we plot the squared RBLM coherent phonon am- 
plitude Q(t) 2 as a function of carrier relaxation time. The 
relaxationless limit is shown as a blue dashed line for com- 
parison. For a relaxation time of 70 fs the squared amplitude 
Q(tf is reduced by 10%. 



carrier relaxation time. The square of Q(t) is shown in 
Fig. [7j We find that for relaxation times greater than 
200 fs, we recover the results obtained in the absence of 
carrier relaxation. This is not surprising since the RBLM 
period is 141 fs and the point about which Q(t) oscil- 
lates is able to follow slow changes in S(t) adiabatically 
while the amplitude of the Q(t) oscillations at long times 
remains unchanged. For relaxation times less than the 
RBLM period, the amplitude of Q(t) at long times de- 
creases, but even for relaxation times slightly less than 
the RBLM period the errors are not too great. For a re- 
laxation time of 70 fs, we see that the square of the Q(t) 
oscillation amplitude and hence the CP intensity is only 
reduced by 10%. Since 70 fs is the shortest photoexcited 
carrier relaxation time measured in the pump-probe ex- 
periments by Dawlaty et al£Z- on graphene samples, we 
feel that our neglect of carrier relaxation effects is a fair 
approximation. Henceforth, we will neglect the effects of 
carrier relaxation for simplicity. 

Since coherent phonon spectroscopy gives direct phase 
information on the coherent phonon amplitude, it is in- 
structive to examine 5 max as a function of pump photon 
energy. This is done in Fig.[5J Fig.JSJa) is the power spec- 
trum of Q(t) at the RBLM frequency. In Fig. [SJ^b) we 
plot S max as a function of pump photon energy. For com- 
parison, the absorption coefficient is plotted in Fig. [5Jc) . 
Near the band edge, we see from Fig.[5Jb) that the pump 
light is strongly absorbed at the V2C1 and V3V1 peaks. 
The resulting increase in the photoexcited carrier den- 
sity increases the coherent phonon driving function and 
enhances the coherent phonon oscillation amplitudes. In 
other words the coherent phonon driving function near 
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FIG. 8: (color online) For 7 zGNR nanoribbon excited by 
Gaussian a laser pulse with polarization vector parallel to 
ribbon length, we plot (a) the coherent phonon power at the 
RBLM frequency (29 meV), (b) the value of S ma x, and (c) the 
initial absorption spectrum as a function of photon energy. 



the band edge is determined by the strength of optical 
absorption between the lowest few hole bands and the 
localized edge states VI and CI. At energies above 3 eV, 
Smax changes sign and the nanoribbon initially contracts. 



D. Detection of coherent phonons 

The generation of coherent phonons in zGNR nanorib- 
bons results in macroscopic oscillations of the carbon 
atoms. The time-dependent macroscopic displacement of 
each carbon atom can be described as a sum of terms each 
of which is proportional to a different coherent phonon 
amplitude^ These displacements in turn modulate the 
optical properties of the nanoribbon through changes in- 
duced in the electronic band structure. 

We illustrate our results in Fig. ® where we simulate 
the results of a hypothetical two-color pump probe ex- 
periment on a 7 zGNR nanoribbon. As usual the fluence 
is 5 x 10 4 J/cm 2 and we assume the pump and probe 
are both linearly polarized with electric polarization vec- 
tor parallel to the nanoribbon length. As indicated in 
Fig. [§ta) the pump photon energy is 2.24 eV which coin- 
cides with the peak in the V3V1 absorption feature. The 
time-resolved differential gain measured by the probe is 
given by 

Ag(fiw,t) = -(a(t}u,t)-a(hLJ,t^-oo)) (10) 
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FIG. 9: (color online) For 7 zGNR nanoribbon excited by 
Gaussian laser pulse with polarization vector parallel to rib- 
bon length, we plot (a) absorption coefficient as function of 
photon energy and (b) differential gain as a function of time 
and photon energy after the transient pump pulse dies out. 
The pump photon energy is taken to be 2.24 eV 



where a{fvjj,t) is the time dependent absorption coeffi- 
cient defined in Eq. ((HJ. The differential gain is shown 
in Fig. [Spa) as a function of probe delay and probe pho- 
ton energy. There are large differential gain oscillations 
near the V2C1 and V3V1 absorption peaks which oscil- 
late at the RBLM period of 0.141 ps, so it is clear that 
the dominant contribution to the time-dependent differ- 
ential gain comes from RBLM induced absorption mod- 
ulation. We find that the ribbon width initially expands 
and begins oscillating about a new and wider equilibrium 
width Wo + AVK where Wo is the initial width before pho- 
toexcitation. When the laser pulse duration is much less 
than the coherent phonon period, the ribbon width os- 
cillates between W and Wo + 2 AW. When the ribbon 
width W(t) = Wo the differential gain vanishes. When 
W(t) — Wo + 2 AW the absorption peaks involving tran- 
sitions to the localized edge states shift to higher energies 
giving rise to positive differential gain on the low energy 
side of the absorption peak and negative differential gain 
on the high energy side. This is clearly seen in Fig.[9jb). 

The experiments we model are degenerate pump/probe 
measurements (i.e. the pump and probe have the same 
wavelength) and the pump and probe excitation en- 
ergy is scanned. After background subtraction, the 
time-dependent differential transmission is Fourier trans- 
formed to obtain the power spectrum. Peaks in this 
power spectrum occur at the coherent phonon frequen- 
cies. Plotting the power spectrum at the RBLM fre- 
quency as a a function of pump/probe photon energy 
gives us the RBLM CP power spectrum. The RBLM CP 
power spectrum is shown in Fig. Efa) as a function of 
pump/probe photon energy. We see strong double peaks 
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FIG. 10: (color online) For 7 zGNR excited by Gaussian 
laser pulse with the pump/probe polarization vector rotated 
an angle from the nanoribbon axis, we plot (a) the RBLM CP 
spectrum as a function of pump/probe energy for several po- 
larization angles from the 50 degrees and (b) the integrated 
RBLM CP intensity for V2C1 and V3V1 transitions. 



in the RBLM CP spectra corresponding to the V2C1 and 
V3V1 transitions. Integrating the area under these dou- 
ble peaks, we obtain the RBLM CP intensity for each 
transition. 

We also investigated the in-plane polarization depen- 
dence of the coherent phonon spectra. Here the polariza- 
tion of the electric field of both the pump and probe pulse 
are rotated by an angle 9 with respect to the nanorib- 
bon length. The dependence of the RBLM CP power 
spectrum on pump/probe polarization angle is shown in 
Fig. HUT a) where the CP power spectra are shown for po- 
larization angles between and 50 degrees. At a polar- 
ization angle of degrees, there is a strong V3V1 signal 
at 2.25 eV and a weaker V2C1 signal at 1.65 eV. As the 
polarization angle increases from degrees to 50 degrees, 
the CP signal becomes weaker. The CP intensity is plot- 
ted on a log scale in Fig. ITDf b) for the V2C1 and V3V1 
transitions. As the polarization angle varies from to 90 
degrees, the CP intensity is strongly quenched. 

It is useful to plot the CP intensity as a function of 
nanoribbon width. We fix the pump and probe elec- 
tric polarization vectors to be parallel to the nanoribbon 
width as this is the polarization for which the CP inten- 
sity is greatest. In Fig. ITTT aL the RBLM CP intensity 
for the V2C1 and V3V1 transitions as a function of the 
number of carbon dimers in the zigzag nanoribbon unit 
cell is plotted against the left axis and the RBLM fre- 
quency is plotted against the right axis. In Fig. Illf b). 
we plot the coherent driving function amplitude S max as 
a function of N a b for the V2C1 and C3V1 transitions. 
We have studied all zGNR nanoribbons for N a b ranging 
from 6 to 17 and find similar results. The driving func- 
tion is positive for low energies (up to just below 3 eV) 
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FIG. 11: (color online) For zigzag nanoribbons excited by 
Gaussian laser pulse with pump and probe polarization vec- 
tor parallel to ribbon length, in (a) we plot the RBLM CP 
intensity for the V2C1 and V3V1 transitions as a function of 
N a b, the number of carbon dimers in the zigzag unit cell, on 
the left axis. On the right axis in we plot the RBLM fre- 
quency ui in eV. In (b) we plot S m ax for the V2C1 and V3V1 
transitions. The ribbon width for each value of N ao can be 
read from the upper axis. 



and then becomes negative for all cases. 



IV. ARMCHAIR NANORIBBON RESULTS 

A. Bandstructure and absorption spectra 

In armchair nanoribbons (aGNR), the nanoribbon 
translational unit cell length is L = y/3 a = 4.31 A and 
the nanoribbon width is W = (N a b — 1) # . Armchair 
nanoribbons belong to one of three families depending on 
the mod numbern = mod( N a b, 3) i 25 ' 35 Following the con- 



vention of Refs. [25| and [35l we label the armchair nanorib- 
bons as cv-aGNR, /3-aGNR, and 7-aGNR families for mod 
2, mod 0, and mod 1 nanoribbons respectively. 

Example bandstructure calculations for each nanorib- 
bon family are shown in Fig. [T2l In our ETB model, all 
armchair nanoribbons have a direct gap at k = 0. From 
the 6 /3-aGNR mod and 7 7-aGNR mod 1 nanoribbon 
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FIG. 12: (color online) Electronic band structure for 
three armchair nanoribbons obtained from the ETB model. 
We have one representative from each mod(iV a (,,3) armchair 
nanoribbon family: (a) 6 /3-aGNR mod semiconducting 
nanoribbon, (b) 7 7-aGNR mod 1 semiconducting nanorib- 
bon, and (c) 8 a-aGNR mod 2 metallic nanoribbon. 



bands in Fig. [T27a) and (b) we see that (3- and 7-aGNR 
ribbons are semiconducting with parabolic bands near 
the band edge. From the 7 cv-aGNR nanoribbon band 
in Fig. [T2Tc) we see that the a-aGNR ribbons have neg- 
ligible band gaps and a linear band edge dispersion re- 
lation. Unlike the case in zigzag nanoribbons, there are 
no localized edge states near the band edge. Armchair 
nanoribbons have direct gaps that arise from quantum 
confinement and edge effects and all the electronic wave- 
functions near the band edge are distributed throughout 
the width of the ribbon. 

In Fig. [13] we plot the direct band gap E g in the ETB 
model as a function of N a b, or equivalently the ribbon 
width W, for the a, /?, and 7 armchair nanoribbon fam- 
ilies. In the case of the a-GNR family, the band gap is 
multiplied by a factor of 10 for clarity. We find that a- 
aGNRs have a small band gap that decreases with ribbon 
width. The j3- and 7-aGNRs have large band gaps that 
decrease as a function of ribbon width, but they fall on 
different curves as seen in the figure. Our ETB results 
are in decent agreement with first principles results re- 
ported in Ref. l25l These results differ from simple tight 
binding (STB) theory which predicts that band gaps in 
a-aGNRs vanish identically and that band gaps for j3- 
and 7-aGNRs as a function of ribbon width fall on the 
same curved. For each family, we fit the band gaps to 
an expression of the form E g = A/(W + B) where W 
is the ribbon width in Angstroms and E g is the band 
gap in eV. For wide ribbons E g A/W. For the range 
6 A < W < 25 A, fitting parameters for each family 
are obtained. For a-aGNRs A = 0.486 eV-A and B = 
3.014 A, for /3-aGNRs A = 13.8 eV-A and B = 2.082 
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FIG. 13: (color online) Direct band gaps for armchair 
nanoribbons in the ETB model as a function of the number 
of dimers N a b in the translational unit cell. Band gaps are 
shown for a-aGNR mod 2, /3-aGNR mod 0, and 7-aGNR mod 
1 families. The band gaps for the a-aGNR mod 2 nanorib- 
bons have been multiplied by a factor of 10 for clarity. The 
nanoribbon width is plotted on the top axis. 



A, and for 7-aGNRs A = 14.97 eV-A and B = 2.286 
A. These results are similar to those obtained in Ref. [35| 
using a semiempirical extended Hiickel theory. The arm- 
chair nanoribbon band gaps in Ref. |35| are described by 
Eg = A/W where values of A are 0.4, 8.6, and 10.4 eV- A 
for a- ft- and 7-aGNRs respectively. 

Since there are three distinct armchair nanoribbon 
families a, ft and 7, we've computed room temperature 
(T = 300 K) absorption spectra for an undoped nanorib- 
bon in each family as a function of photon energy for 
light with electric polarization vectors parallel to the rib- 
bon length. In computing the absorption coefficient us- 
ing Fermi's golden rule (Eq. (|8|)) we assume Lorentzian 
lineshapes with the same 0.2 eV FWHM used in Fig. [3] 

In Fig. ll4f a) we plot the absorption spectrum in an un- 
doped 6 /3-aGNR mod semiconducting nanoribbon for 
light polarized parallel to the nanoribbon length. The 
total absorption spectrum, shown as a solid black line, is 
a sum of absorption spectra due to three major interband 
direct transitions labeled En, E22 and -E33. These transi- 
tions are shown in the band diagram in the right panel of 
Fig. lT4T a). The sharp peaks in the spectra are due to von 
Hove singularities at each of the transition band edges. 
The room temperature Fermi level is Ep = —29.5 meV. 
The Fermi level is shifted slightly to the valence band 
edge since the conduction and valence bands are asym- 
metric in the ETB model. Sasaki et al£2- have studied 



optical selection rules in armchair nanoribbons and find 
that for light polarized parallel to the nanoribbon length, 
the allowed interband transitions are direct transitions 
between valence band i and conduction band i. In order 
to conform to a widely used convention in the literature, 
the transitions in Fig. IT47 a) are labeled En as opposed 
to the ViCi notation used in Section IIII1 We see that the 
optical selection rules in armchair nanoribbons are quali- 
tatively different from those in zigzag nanoribbons where 
direct transitions are forbidden for polarization parallel 
to the ribboni22i22 

In Fig. IT47 b) we show similar results for an undoped 
7 7-aGNR mod 1 semiconducting nanoribbon. Again 
the total absorption spectrum is the sum of absorption 
spectra due to three interband direct transitions En, E22 
and E33 and the room temperature Fermi level Ep = 
—32.7 meV is shifted slightly towards the valence band 
edge. 

In Fig. [T4T c) we show results for an undoped 8 a-aGNR 
mod 2 metallic nanoribbon. In the left panel of Fig. IT47 c) 
we plot the absorption spectrum for light polarized paral- 
lel to the nanoribbon length. The total absorption spec- 
trum, shown as a solid black line, is a sum of absorp- 
tion spectra due to three major interband direct transi- 
tions labeled En, E22 and E33. In the metallic armchair 
nanoribbons, we follow convention and label optical tran- 
sitions starting with the lowest one defined as En- For 
clarity, transition energies are sometimes labeled or 
E™ in semiconducting and metallic nanoribbons, respec- 
tively. Our notation should cause no confusion since it 
is clear from the context whether we are talking about 
semiconducting or metallic nanoribbons. From the band 
diagram in Fig. [T4Tc) we see that En corresponds to a 
direct transition between the second electron and hole 
bands. Optical transitions between the first electron and 
hole bands with linear dispersion don't contribute to the 
absorption spectrum. They do, however, pin the Fermi 
level as a function of T at E-p = meV as can be seen in 
Fig. [He). 



B. Phonon dispersion relations 

Phonon dispersion relations for 6 /3-aGNR mod 0, 
7 7-aGNR mod 1, and 8 a-aGNR mod 2 nanoribbons 
obtained in our VFF model are shown in Figs. H~5T a). 
(b) and (c), respectively. As in the zigzag nanoribbon 
case, dispersion relations and normal mode displacements 
are obtained by solving the dynamical matrix eigenvalue 
problem. The in-plane and out-of-plane modes can again 
be treated independently. There are 2 N a b out-of-plane 
modes (red dash-dotted lines) and 4 N a b in-plane modes 
(black solid lines). In all cases, we find that N a b of the 
in-plane q — modes can be coherently excited. The 
lowest of these modes is an RBLM mode. The RBLM 
phonon dispersion relations are shown in Fig.[TS]as thick 
blue lines. 

The q = RBLM phonon energies in aGNR nanorib- 
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FIG. 14: (color online) Band absorptions with electric polarization vector parallel to the ribbon and band diagram showing 
the room temperature Fermi level Ep with transitions involved in the three lowest absorption peaks for (a) 6 /3-aGNR mod 
nanoribbon, (b) 7 7-aGNR mod 1 nanoribbon, and (c) 8 a-aGNR mod 2 nanoribbon. Total absorption is the sum of absorption 
due to several transitions. In the band diagrams, blue lines are valence bands and red lines are conduction bands. 




q (nil) q (n/L) q (n/L) 



FIG. 15: (color online) Phonon dispersion relations for 
aGNR nanoribbons calculated in the valence force field model 
for (a) 6 /3-aGNR mod 0, (b) 7 7-aGNR mod 1, and (c) 
8 a-aGNR mod 2 nanoribbons. The out-of-plane modes are 
red dash-dotted lines while the remainder are in-plane modes. 
The dispersion band associated with the q — RBLM mode 
is shown as a thick blue line. 



bons obtained from our VFF model can be fit to fkc = 
AW P + B where W is the armchir nanoribbon width in 
Angstroms and hio is the q = RBLM phonon energy in 
eV. The fitting parameters A = 0.31376, P = -0.95579 
and B = 0.00328 are obtained from a least squares fit 
to our VFF results for 4 < N a f, < 25 corresponding to 
widths in the range 4A<W<30A. 



C. Generation of coherent phonons 

In our simulation, we excite coherent RBLM phonons 
with a 20 fs Gaussian laser pulse. The fluence F is 



5 x 10 4 J/cm and the pump and probe pulse are lin- 
early polarized with an electric polarization vector par- 
allel to the nanoribbon length. Time dependent distri- 
bution functions are obtained by integrating the photo- 
generation rates in Eq. © using the equilibrium dis- 
tribution functions as initial conditions. We assume a 
FWHM pump spectral width of 50 meV which appears 
in the broadened delta function in Eq. ([3]). Using the 
time-dependent distribution functions, the RBLM coher- 
ent phonon driving function S(t) can be obtained from 
Eq. ©. 

We examine 5 max as a function of pump photon en- 
ergy and our results for 6 /3-aGNR mod semiconduct- 
ing nanoribbons are shown in Fig. I16f a) where S ma-X is 
shown as a function of pump photon energy. For com- 
parison, the absorption coefficient is also plotted in the 
lower panel of Fig. ITBT a) . 

Near the band edge, we see from Fig. ITBT a) that the 
pump light is strongly absorbed near the En and E22 
peaks. The resulting increase in the photoexcited car- 
rier density increases the coherent phonon driving func- 
tion and enhances the coherent phonon oscillation ampli- 
tudes. Photoexcitation by the pump causes the nanorib- 
bon to initially expand for pump photon energies near the 
En transition and to initially contract for pump photon 
energies near the E22 transition. We find this to be true 
for all /3-aGNR mod semiconducting nanoribbons. 

Qualitatively different results are obtained for 7-aGNR 
mod 1 nanoribbons. In Fig. fTBTb) we plot S max as a 
function of pump photon energy for a 7 7-aGNR mod 1 
nanoribbon and find that photoexcitation by the pump 
causes the nanoribbon to initially contract for photon en- 
ergies near the En peak and initially expand for photon 
energies near the E22 peak. This is found to be true for 
all 7-aGNR mod 1 semiconducting nanoribbons. 

In Fig. fTBTc) we show results for an 8 a-aGNR mod 
2 metallic nanoribbon excited by a laser pulse polarized 
parallel to the ribbon length. From Fig. ITBT c). we see 
that photoexcitation by the pump causes the nanorib- 
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FIG. 16: (color online) The coherent phonon power, the value of Smax, and the initial absorption spectrum are plotted as a 
function of photon energy for (a) 6 /3-aGNR mod nanoribbon (RBLM frequency = 59 meV) , (b) 7 7-aGNR mod 1 nanoribbon 
(RBLM frequency = 51 meV), and (c) 8 a-aGNR mod 2 nanoribbon (RBLM frequency = 44 meV). The excitation is due to a 
Gaussian laser pulse with pump and probe polarization vectors parallel to the ribbon length 




Polarization Angle (degrees) 



polarization angle making an angle 9 with the nanorib- 
bon axis. Our results are shown in Fig. [T7] and we can 
see that CP intensity for the En and E22 features in 
7 7-aGNR mod 1 semiconducting nanoribbons drop off 
sharply as a function of polarization angle. 

We simulated CP intensity as a function of N a b in arm- 
chair nanoribbons photoexcited by ultrafast Gaussian 
laser pulses polarized parallel to the nanoribbon length. 
The results are shown in Fig. [PBT a) for the En and E22 
features in a-, j3-, and 7-aGNR nanoribbons. The corre- 
sponding values of the coherent phonon driving functions 
Smax are shown in Fig. [PBT b). 



V. ANALYSIS OF THE ELECTRON-PHONON 
INTERACTION 



FIG. 17: (color online) Integrated RBLM CP intensity for 
En and E22 transitions excited by a Gaussian laser pulse as 
a function of pump/probe polarization angle in a 7 7-aGNR 
mod 1 nanoribbon. 



bon to initially expand for photon energies near the En 
transition. For photon energies near the E22 transition, 
the situation is more ambiguous. 

We have studied the CP intensity for ultrafast photoex- 
citation near the En and E22 transitions in a 7 7-aGNR 
mod 1 semiconducting nanoribbon as a function of polar- 
ization angle. In our simulation, we assume the polariza- 
tion of pump and probe are parallel with the pump/probe 



Following our previous study, 9 in this section we exam- 
ine the fc-dependent electron-phonon interaction in the 
effective mass approximation to explain why some GNRs 
start their coherent RBLM oscillations by initially ex- 
panding while others start the oscillations by initially 
shrinking. In the present discussion we will focus our 
attention on the aGNRs, in which we can directly use 
the wavefunctions formulated in previous papersi^iSi In 
the case of zGNRs, we have to consider a special local- 
ized wavefunction due to the presence of edge states at 
which the En transition occurs. 37 Using such a wave- 
function, we obtain a constant electron-phonon matrix 
element that does not depend on mod (N a i,, 3) of the 
zGNRs, and thus is consistent with our results in Sec. lIIII 
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FIG. 18: (color online) For a-, /?-, and 7-aGNR nanoribbons 
excited by Gaussian laser pulses with polarization vector par- 
allel to ribbon length, we plot (a) RBLM CP intensity for 
Eu and E22 transitions and (b) the corresponding coherent 
phonon driving functions S ma x as functions of N a b, the num- 
ber of carbon dimers in the zigzag unit cell. The curves with 
upward pointing arrows are for the Eu transitions and the 
downward pointing arrows are for the £22 transitions. The 
ribbon width for each value of N a t can be read from the upper 
axis. 



The details will be presented elsewhere. 

The electron-phonon matrix element M cp is a sum 
of conduction band (c) and valence band (v) electron- 
phonon matrix elements, which represent the electron 
and hole contributions, respectively, 



M ep = Mg p 



M, 



cp 



= (c|J3ep|c) ~ (V|if e: 



(11) 



where H ep is the GNR electron-phonon interaction 
Hamiltonian and the minus sign comes from the opposite 
charges of electrons and holes. In a nearest-neighbor ef- 
fective mass approximation, the RBLM H cp for an aGNR 
can be written as 



H cp = u. 



a —828. 

yon 2 
arm | g fi 

9 yon 



(12) 



where g on (g g) is the on-site (off-site) coupling constant 
in eV, while it arm is a ribbon width- or A^-dependent 
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FIG. 19: Cutting lines for (a) mod /3-aGNR and (b) mod 1 
7-aGNRs near the Dirac K point. To make clear the definition 
of O(k), in this figure 0(k) is shown for an arbitrary k at 
Eu. In fact, in the case of mod and mod 1 aGNRs the 61 
transitions occur at O(k) = and O(k) = n, respectively. 
The difference between the mod and mod 1 aGNRs can be 
understood from the position of the Eu or E22 cutting lines 
relative to the K point. 



phonon amplitude. See Appendix [Dl for the derivation of 
this Hamiltonian. Next, to obtain M ep in Eq. ([TTj) . we 
use the following wavefunctions 
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~.ik-r 



e -i0(k)/2 

c +ie(k)/2 



■ikr 



e -iS(k)/2 

_ p +ie(k)/2 



for conduction and valence states, respectively. In 
Eq. (|13p . S is the surface area of graphene and 8(k) is an 
angle of k = (k x , k y ) measured from the fc^-axis (ribbon 
width direction). These wavefunctions are suitable near 
the graphene Dirac K point and thus they can explain 
well the aGNR lattice response especially at relatively 
low energy En and E22 optical transitions. 

By inserting the wavefunctions in Eq. (|13[) into 
Eq. pT|) . we obtain 



(c|iT ep |c 



MH, 



cp 



Warm 
v) = U arm 



(2y-c.ii - ffoff cos 0(k)) , (14a) 
(2g n + <7offcose(k)), (14b) 



and thus 



M, 



"arm (-2g S C0s6(k)) 



(15) 



From this equation, we can analyze the N a b and En de- 
pendence of the aGNR initial lattice response. First of 
all, we should note that g g (which is usually taken to be 
6.4 eV as mentioned in Ref. H) and u arm (cf. Eq. (|D6p ) 
are always positive, while cosO(k) can either be positive 
or negative depending on the value of k at which the Eu 
transition occurs. 

According to the definition of the driving force in 
Eqs. ([I} and ^j, a negative (positive) M cp value cor- 
responds to a positive (negative) 5 max - Therefore, a 
positive (negative) cos0(k) is related to a contraction 
(expansion) of the ribbon width. Using this argument, 
we can classify the aGNR lattice response based on the 
aGNR types. For example, let us consider semiconduct- 
ing mod /3-aGNR and mod 1 7-aGNRs. The cutting 
line position for their En and E22 optical transitions are 
just opposite to each other. For a mod /3-aGNR, we see 
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FIG. 20: (color online) Driving force S max and initial ab- 
sorption spectrum as a function of photon energy for several 
mod semiconducting /3-aGNRs: (a) N a b = 6, (b) N a b = 9, 
(c) N a b = 12, and (d) N ab = 15. Positive (negative) 5 max 
at En (E22) corresponds to an expansion (contraction) of 
the ribbon width. 8E — E22 — En is also shown, which is 
decreasing as a function of N a b- 



FIG. 21: (color online) Driving force S max and initial ab- 
sorption spectrum as a function of photon energy for several 
mod 1 semiconducting 7-aGNRs: (a) N a b = 7, (b) N a b = 10, 
(c) N ab = 13, and (d) N ab = 16. Negative (positive) 5" max 
at £n (E22) corresponds to a contraction (expansion) of the 
ribbon width. 



that cos0(k) becomes positive (negative) at En (E22), 
and thus the aGNR starts the coherent phonon oscilla- 
tions by expanding (shrinking) its width. This can be 
seen in the illustration of 9(k) in Fig. [T5] The opposite 
behavior is true for mod 1 7-aGNRs. These arguments 
are consistent with our numerical results in Figs. [2"01 and 
[2T1 in which we show the aGNR driving force trends 
within the mod (3 family and mod 1 7 family, respec- 
tively. 

However, the driving force trends for the mod 2 metal- 
lic a-aGNRs (see Fig. |2"2"|) cannot be explained nicely by 
the effective mass theory for several reasons. The main 
reason is that, in the metallic a-aGNRs, there are two 
cutting lines with the same distance from the K point, 
which can be assigned as the lower and higher branches 
of an En transition. Both branches contribute to a spe- 
cific En and thus we have to sum up the matrix elements 
from each contribution to obtain M cp . For example, if 
the ID fc-points for the lower and higher branches of En 
are the same, the matrix elements will cancel each other 
because cos 8(k)+cos(7r— O(k)) = 0. In this case, the CP 
amplitude will be generally small for the mod 2 metallic 
aGNRs compared to the mod or mod 1 semiconducting 
aGNRs. In the real case, we always have slightly differ- 
ent /c-points for the lower and higher branches due to the 
trigonal warping effect^ from which the nonzero M ep 
value gives information about an expansion or contrac- 
tion of the ribbon width. 

Near the En transition, the metallic cv-aGNR initial 
lattice response is always an expansion for all N ab . On 
the other hand, near the E22 transition, the response is 
expected to be always a contraction, though we see in 
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FIG. 22: (color online) Driving force S ma x and initial absorp- 
tion spectrum as a function of photon energy for several mod 2 
metallic aGNRs: (a) N ab = 8, (b) N ab = 11, (c) N ab = 14, 
and (d) N ab = 17. Positive (negative) S max at En (E22) cor- 
responds to an expansion (contraction) of the ribbon width. 



Figs. [22] (b)-(d) the trends is not true for larger N a b. We 
notice that the difference Se — E22—E11 might determine 
whether or not the lattice response at the E22 feature 
of a given mod 2 N a f, a-aGNRs will clearly follow our 
effective mass theory. We guess that if Se is large enough 
(ps 2eV as in the 8 a-aGNR), the lattice response at E22 
should not be ambiguous. However, Se decreases with 
increasing N a b as can be understood from a cutting line 
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TABLE I: Initial lattice behavior due to coherent phonon os- 
cillations at _En and E22 in aGNRs. 



family 


En 


E22 


mod 


expand 


contract 


mod 1 


contract 


expand 


mod 2 


expand 


expand or contract 



argument. Thus the lattice response at E 22 for the 
larger mod 2 aGNRs becomes opposite to that for the 
smaller mod 2 aGNRs (e.g. N ab = 2, 5, 8). 

We should also note that the En and E 2 2 values of 
the metallic a-aGNRs are close to the order of £33 and 
E44 values for the semiconducting aGNRs. This gives 
another reason why the effective mass theory cannot ex- 
plain the metallic a-aGNR trends. In this energy region, 
the nearest-neighbor effective mass theory should be ex- 
tended to include longer-range nearest-neighbor interac- 
tions, 3 ' 1 which is beyond the scope of this paper. Never- 
theless, the present discussion on the effective mass the- 
ory has already given physical insight for explaining the 
initial lattice response of mod and mod 1 semiconduct- 
ing aGNRs. Even in the case of mod 2 metallic a-aGNRs 
we can see that the effective mass theory is still able to 
explain the lattice behavior near the En region. Finally, 
we summarize the lattice behavior at En and E 22 tran- 
sitions for all families of aGNRs in Table HI 



VI. SUMMARY 

We have developed a microscopic theory for the gen- 
eration and detection of coherent phonons in graphene 
nanoribbons using an extended tight-binding model for 
the electronic states and a valence force field model for 
the phonons. We consider zigzag and armchair ribbons 
denoted N a b zGNR and N a b aGNR, respectively, where 
N a b is the number of AB carbon dimers in the transla- 
tional unit cell. In coherent phonon spectroscopy, ultra- 
fast laser pulses generate electrons and holes in the con- 
duction and valence bands of a graphene nanoribbon. If 
the pulse duration is less than the phonon oscillation pe- 
riod, the photogenerated carriers couple to the phonons 
through the deformation potential electron-phonon in- 
teraction and the nanoribbon lattice undergoes coherent 
macroscopic phonon oscillations. The coherent phonon 
amplitudes satisfy a driven oscillator equation (derived 
from the Heisenberg equation) with a driving function 
that depends on the electron-phonon interaction matrix 
elements and the photoexcited carrier distribution func- 
tions. For large laser spot sizes, only q — phonon modes 
can be coherently excited and these modes are said to be 
CP active. Coherent phonons are detected using a de- 
layed probe pulse to measure the time dependent oscilla- 
tions in the differential transmission. Taking the Fourier 
transform of the differential transmission with respect to 
probe delay time, we obtain the coherent phonon spec- 
trum as a function of phonon frequency with peaks in the 



CP spectrum corresponding to excited coherent phonon 
modes. 

For both zGNRs and aGNRs there are A^;, CP active 
modes that vibrate in the plane of the nanoribbon. In 
all cases, the CP active mode with the lowest frequency 
is the RBLM mode. For 20 fs Gaussian laser pulses we 
ignore carrier relaxation effects and integrate Eq. (J3)), the 
photogeneration rate obtained from Fermi's golden rule, 
to obtain the photoexcited carrier distribution functions. 
Using these carrier distribution functions we obtain the 
driving function for RBLM coherent phonons as a func- 
tion of the pump photon energy and polarization angle. 
In all cases we find that the driving function .SVax imme- 
diately after photoexcitation is greatest for light polar- 
ized parallel to the nanoribbon axis. For photoexcitation 
near the optical absorption edge, we find that the coher- 
ent phonon driving function for the RBLM mode is much 
larger for zigzag nanoribbons where An ^= transitions 
involving localized edge states provide the dominant con- 
tribution to the CP driving function. The sign of •S'max 
is interesting since it gives phase information that can be 
measured in CP spectroscopy. In zigzag nanoribbons, the 
ribbons initially expand for pump photon energies near 
the V2C1 and V3V1 absorption features. In armchair 
nanoribbons the phase behavior is different for the a, /3 
and 7 families where An = transitions between zone 
folded quantum confined states give rise to the CP driv- 
ing function. In /3-aGNR mod semiconducting aGNRs 
the ribbon initially expands for pumping near the En 
absorption peak and contracts for pumping near E 22 ■ In 
7-aGNR mod 1 semiconducting aGNRs the ribbon ini- 
tially contracts for pumping near the En absorption peak 
and expands for pumping near E 22 . In the case of the 
metallic a-aGNR mod 2 ribbons, the ribbon is seen to 
expand for pumping near the En feature while the be- 
havior near the E 22 feature is ambiguous. Such lattice 
responses are also discussed by considering the electron- 
phonon interaction within an effective mass theory. 
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Appendix A: Electronic states in extended tight 
binding model 

In treating electronic states in zigzag and armchair 
nanoribbons within the ETB model, we assume that all 
bond lengths and bond angles between A and B carbon 
atoms are the same as that in planar graphene. In addi- 
tion, we assume the ribbon edges are unpassivated. Since 
the nanoribbons are assumed to be perfectly planar, there 
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is no mixing of a and it orbitals and the two types of 
bonds can be considered separately. The a bonds ac- 
count for the mechanical properties of the ribbon while 
the 7r bonds determine the optical properties near the 
band edge and play the dominant role in the generation 
of coherent phonons. We treat the valence tt and conduc- 
tion 7r* bands using the third-nearest-neighbor extended 
tight binding (ETB) model developed by Porezag et al. 
for carbon compounds^ The Hamiltonian and overlap 
matrix elements between tt orbitals on different sites are 
parameterized by analytic expressions that depend on 
the C-C bond length. These parameterizations were cal- 
ibrated using DFT results for a local orbital basis set in 
the local density approximation (LDA) for a wide range 
of carbon compounds. 

A graphene nanoribbon only has translational sym- 
metry along its length with period L so the Brillouin 
zone is one dimensional with a wavevector k satisfying 
\k\ < n/L. The tight binding wavef unction for a tt elec- 
tron with wavevector k is given by 

B N ab oo 

M*) = J2 E<^-( fc ) E eWM*-K) ( A1 ) 



s=A j=l 



l= — oo 



where (f>Tv(r) are atomic 2p z orbitals, R SJ - is the position 
of carbon atom sj in the 1-th translational unit cell. 

The tight binding expansion coefficients C S j{k) in 
Eq. (jAll) satisfy the generalized eigenvalue equation 

£ H$ (k) C slj , (k) = E{k) S$ (k) C s , y (k). (A2) 



s'j' 



s'j' 



The Hamiltonian matrix elements are given by 



ff#(*) = Ee' (M), *£?( 



R 



S'J' 



R 



S'J' 



(A3) 



where H^^ K {d) is the Hamiltonian matrix element be- 
tween two carbon 2p z tt orbitals as a function of the dis- 
tance d between them. A similar expression holds for 
S?j, (k) with HpJ^(d) replaced by the overlap matrix el- 
ement Sp^(d). The ETB matrix elements vanish at a 
finite cutoff distance, and the sum over I in Eq. (|A3[) 
includes all contributions out to third nearest neighbor 
shells. Analytic expressions for and can be 

found in Ref. d 

Solving Eq. (|A2[) we obtain the electronic states E n (k) 
and wavefunctions where n = 1 . . . 2N a b is the subband 
index in order of increasing energy. The first N a b levels 
are the valence bands and the remainder are conduction 
bands. 



mode h 20 ' 21 ' 40 We let U z SJ (i) be the displacement from 



equilibrium of a carbon atom at R SJ -. 
motion are given by 



The equations of 



dt 



^y^j K-sjl;s'j'l' 
s'j'l' 



U' (i)-U'.(i) . (Bl) 



The force constant tensor K. s jl;s'j'l' depends on the dis- 
tance between atoms sjl the s'j'l'. In order to describe 
bond twisting interactions involving four atoms it is nec- 
essary to retain up to fourth nearest neighbor shell force 
constant tensors. If 77 labels nearest neighbor shells and 
d v is the distance between them, the force constant ten- 
sor connecting 77-th neighbor atoms is 



K 



sjl;s'j'l' 



<t>?)T x Ty 




(<& ~ <t> n r)T: 



r 'y 



x'y 







(B2) 

where effy is the 77-th neighbor force constant parameter 
in the radial bond-stretching direction and <\P ti and <p^ 
are 77-th neighbor force constants for the transverse in- 
plane and out-of-plane directions. For neighbor shells 
77 = 1 ... 4, we use the force constants tabulated in Table I 
of Ref. SU and reproduced in Table 9.1 of Ref. H The 
vector r is a vector of unit length from atom sjl to s'j'l'. 

Taking the Fourier transform of the atomic displace- 
ments, we obtain 



ui,(t) = 



1 



E U ^) 



(B3) 



where XJ S j (q) are the normal mode displacement vectors 
and q is parallel to the nanoribbon axis. The phonon 
energies and mode displacement vectors are obtained by 
solving the dynamical matrix eigenvalue problem 

E V'WiQ) U ^'(9) = Muj{ q f V SJ (q), (B4) 



where M is the carbon atom mass and u(q) is the eigen- 
frequency of the phonon mode. In terms of the force 
constant tensor, the dynamical matrix is given by 

T>s 3 ;s'j'(q) = 5 ss '5jji K s jO;s"j"l" 
s"j"l" 

- E'<— e-^-^.CBS) 
l" 

For each value of q the dynamical matrix eigenvalue 
problem can be solved to obtain phonon energies hu} m (q) 
and displacement vectors U™(g) where m labels the 
phonon modes. 



Appendix B: Phonon modes in valence force field 
model 

We treat lattice dynamics in zigzag and armchair 
nanoribbons using a Born- von Karman valence force field 



Appendix C: Electron-phonon interaction in ETB 
model 

The second-quantized electron-phonon interaction in 
zigzag and armchair nanoribbons in the ETB model is 
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obtained by evaluating the integral 

H ep = J dr $t(r) H ep (r) #(r) (CI) 

where the deformation potential electron-phonon inter- 
action Hamiltonian is 

sjl 

The electron field operator ip(r) is 

VK r ) = ^c nk ip nk (r), (C3) 

where ipnk(?) are the tight binding wavefunctions for 
wavevector fc in Eq. (|A1[) with n = 1 . . . 2N a b labeling 
the subband levels. The annihilation operator for this 
electronic state is c nk . 

After carrying out this integral over the electronic co- 
ordinate, we obtain the final result by replacing the classi- 
cal atomic displacements with the second-quantized 
operator 

(C4) 

Here p is the mass density per unit length and Lq is the 
nanoribbon length. The atomic displacements e™ (q) are 
eigenvectors of the dynamical matrix in Eq. (|B5|i . They 
satisfy the normalization condition 

5>3(9)*' e 3' (9) (C5) 
si 

The final result for the second-quantized electron- 
phonon interaction is 

Hep = Y^ E M ^k';nk C l'k' C nk (bmq + b^g), (C6) 
m n'k';nk 

where q = k' — k. The interaction matrix elements is 
given by 

M%l, ;nk = -A m (<l) E C:, r (n'k')C S3 (nk)x 

sj;s'j' 

E e lL(fc '- fc ' r) e?,,.,, (g) • X(s'j'l'; s"j"; sjl),(C7) 

U';s"j" 



where A m (q) — hj \/2pL^ huj m (q) is the quantized 
phonon amplitude. In Eq. (|C7j) L is the length of the 
translational unit cell, C s j(nk) are the electron tight- 
binding expansion coefficients in Eq. (|A1|) , e™ (q) are the 
atomic normal mode displacements obtained by diago- 
nalizing the dynamical matrix in Eq. (|B5[) . and A is the 
deformation potential vector. 



The deformation potential vector is the three center 
integral 

\ = J dr^Cr-R^OV^r-R^O^Cr-Ri,.)- (C8) 

We evaluate A using the 2p z atomic wavefunctions and 
screened atomic potential for carbon in Ref. [22| obtained 
from an ab initio calculation in graphene^ The ir or- 
bitals are expanded as 

<^(r) =z Y,h ex P (-^2 ) . ( C9 ) 
1=1 ^ ° l ' 

where the z direction is perpendicular to the nanorib- 
bon plane. Similarly the screened atomic potentials are 
expanded as 




(CIO) 



Values of vi, 07 and 77 are tabulated in Table I of 
Ref. [22]. Substituting the expansions (|C9I) and (|C10|) 
into Eq. (|(J8[) we obtain an expansion for A in terms 
of three-dimensional integrals which can be evaluated 
analyticall y 22 ' 23 



Appendix D: RBLM electron-phonon interaction in 
effective mass model 

In this section, we derive Eq. (fT2]) which gives the 
RBLM electron-phonon interaction Hamiltonian H cp for 
an armchair nanoribbon. The Hamiltonian can be de- 
composed into the on-site and off-site Hamiltonians. The 
details of the on-site and off-site interactions are given in 
Sasaki's work on the deformation-induced gauge field in 
graphene^ We will directly use his results in formulat- 
ing the on-site and off-site Hamiltonians. The on-site and 
off-site interactions are induced by a lattice deformation 
which gives rise to a change in the transfer integral and 
a change in the potential between A and B atoms. We 
adopt a coordinate system shown in Fig. [23] to derive 
H cp . 



On-site Hamiltonian 

The on-site Hamiltonian can be expressed in terms of 
the divergence of ua and ub, which represent the dis- 
placement vector of A-atom and B-atom in the graphene 
unit cell, respectively. This Hamiltonian is written as 

n on = g on (^ " B(r) v .° A ( r ))- ( D1 ) 
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FIG. 23: Upper panel shows displacements of B-atoms at 
r; + Ra (a — 1, 2, 3), that is ub(i"i + R a ), which give rise to 
a deformation potential at A-atom of r;. Lower panel shows 
local modulations of the hopping integral defined by S"/q(t) 
(a — 1,2,3). In this coordinate system we have the nearest- 
neighbor vectors Ri = (0,a cc ), R2 = (— V3/2, — l/2)a cc , 
R 3 = (V^/2,-l/2)occ, 

where a cc = a/-\/3. Here £ — 3a cc /2 

is used in Eq. ()D9[I . 



substituting u(r) in Eq. (|D4[) into Eq. (|D3|) . we obtain 
H oa = . arm (°~ g l) , (D5) 

where, for simplicity, here we have defined 

^arm(^) = km ^rn cos(fc m a;). (D6) 

Note that M ar m is always positive because the atomic po- 
sition x is within the range — W/2 < x < W/2. 

Off-site Hamiltonian 

To derive the off-site interaction Hamiltonian, we 
start with the fact that the lattice deformation modifies 
the nearest-neighbor hopping integral locally as —70 — > 
—70 + £70 (r,-) (a = 1, 2, 3). The corresponding perturba- 
tion of the lattice deformation is given by 



H d 



cform 



E 

ieA a=l,2,3 



E <57o Q (r i )[(cf + J t cf + (cf)t c f +a ] 



(D7) 

where cf" is the annihilation operator for a n electron 
on an A-atom at position 17, and (cS_ a )' is the creation 
operator for a n electron on a B-atom at position 17 +Q (= 
17+R Q ). This perturbation gives rise to scattering within 
a region near the i^T-point of graphene whose interaction 
is given by a deformation- induced gauge field A(r) = 
(A x (r), A y (r)) as vf<J • [p + A(r)], where vf is the Fermi 
velocity, p = — iftV is the momentum operator, and a is 
the Pauli matrix. The deformation-induced gauge field 
A(r) is defined from £7$ (r) (a — 1, 2, 3) as2& 



For the discussion of the RBLM electron-phonon inter- 
action, we rewrite Eq. (jDlj) as follows: 



He. 



/ u A (r) + u B (r) 

/ u A (r) - u B (r) 
V 2 



(D2) 



where g on denotes the gradient of the atomic potential at 
r, (To is the identity matrix, and a z is the z-component of 
the vector of Pauli matrices. For the RBLM in armchair 
nanoribbons we have u A (r) = ub(i") = u(r). Therefore, 
Eq. (|D2|) can be simplified to be 



(D3) 



Hon = 5onCT V ■ u(r). 

The RBLM oscillation can be expressed by 
u(r) = u m sin(k m x)e x , 



(D4) 



where u m is the maximum amplitude of the RBLM oscil- 
lation at the armchair edge. The wavevector k m is then 
specified by ir/W, where W is the ribbon width. By 



1*4,(1-) = ^(r) - - [*yg(r) + 5 7o 3 (r)] , 
t, F A y (r)^^[<57„ 2 (r)-<5 7 o 3 (r)]. 



(D8) 



Similar to the RBM electron-phonon interaction in 
carbon nanotubes^ the perturbation to the nearest- 
neighbor hoppping integral for the RBLM electron 
phonon interaction in the armchair nanoribbons is given 
by 



^(r)^R a .{u(r + R a )-u(r)}, 

ta c .c. 



(D9) 



where g g is the off-site coupling constant and £ = 3a cc /2 
(see the lower panel of Fig. [23)) . Here the displacement 
vector of a carbon atom at r in general is expressed by 
u(r) = [u x (r),u y (r)}. Using a Taylor expansion, we ap- 
proximate Eq. (|D9[) as 



^° Q(r) = if Ra ' {(Ra ' v)u(r)} - 



(DIO) 

Using Ri, R2, and R3 in Fig. [531 we obtain the 
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deformation-induced gauge field of Eq. (|D8I) as follows: be written as 



v F A x (r) 
v F A y {r) 



gog 

2 
5off 



5w x (r) 9%(r) 



9a; dy 
du x {r) du y (r) 



dy 



dx 



(Dll) 



Inserting u(r) in Eq. (|D4[) into this last equation, we 
then obtain 

V F A X = -^-k m U m COs{k m x) = -ytlarm, (D12) 

while v F Ay — 0. Therefore, the off-site Hamiltonian can 



%off = ct x vfA x = u. 



»> t: ! gog Q 



(D13) 



Finally, we can get the electron-phonon interaction 
Hamiltonian of Eq. (fl~2|): 



(D14) 



-ffcp — Hon + T~LoS 



arm I _ 3off 



J/on 2 
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